function num = hwozd_Snms(n,m,s)
    a = (n + abs(m)) / 2;
    b = (n - abs(m)) / 2;
    num	= ((-1)^s)*lans_fac(n-s);
	den	= lans_fac(s)*lans_fac(a-s)*lans_fac(b-s);
    num = num / den;
